A New Cosmological Model of Quintessence and Dark Matter 
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We propose a new class of quintessence models in which late times oscillations of a scalar field give 
rise to an effective equation of state which can be negative and hence drive the observed acceleration 
of the universe. Our ansatz provides a unified picture of quintessence and a new form of dark matter 
we call Frustrated Cold Dark Matter (FCDM). FCDM inhibits gravitational clustering on small scales 
and could provide a natural resolution to the core density problem for disc galaxy halos. Since the 
quintessence field rolls towards a small value, constraints on slow-roll quintessence models are safely 
circumvented in our model. 



PACS number(s): 04.62. +v, 98.80.Cq 

The recent discovery that type la high redshift super- 
novae are fainter than they would be in an Einstein-de 
Sitter universe suggests that the universe may be acceler- 
ating, fuelled perhaps by a cosmological constant or some 
other field possessing long range 'repulsive' effects [|l],p|. 
The acceleration of the universe is related to the equation 
of state of matter through the Einstein equation 
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Pc + Px{i + 3wx)] 



(1) 



for cold matter pc and X-matter with equation of state 
Px — wpx- Clearly a necessary (but not sufficient) con- 
dition for the universe to accelerate is wx < —1/3. In 
other words the equation of state of X-matter must vi- 
olate the strong energy condition (SEC) so that px + 
2>Px < 0. Investigations of cosmological models with 
Om -f Q.X — 1 have demonstrated that these models out- 
perform most others in predicting the correct form for 
the large scale clustering spectrum, accounting for CMB 
anisotropics on large and intermediate angular scales 
and providing excellent agreement with the luminosity- 
distance relation obtained from observations of high red- 
shift supernovae pj. In addition, flat models are com- 
pelling from a theoretical viewpoint since they agree with 
generic predictions made by the inflationary scenario. 

The literature describing phenomenological forms of 
matter violating the SEC is vast (see [Q for a recent re- 
view). Nevertheless two kinds of matter have been sin- 
gled out in recent times as being of special interest: 



1. A cosmological constant Px 

A EE pa'/SttG. 



-px [wx = -1), 
2. A scalar field rolling down a potential V{4)). 



For fields rolling sufficiently slowly </) ~ and Tit ~ 
V{(t))gik, so that ^(0) plays the role of a time-dependent 
A-term. Although appealing, models with the simplest 
potentials including V oc m?(jp' run into problems similar 
to those encountered by a cosmological constant. The 
enormous overdamping of the scalar field equation dur- 
ing radiation and matter dominated epochs causes V{4>) 



to remain unchanged virtually from the Planck epoch 
Zpi ^ 10^® to z ~ 2 H resulting in an enormous differ- 
ence in the scalar field energy density and that of mat- 
ter/radiation at early times. This leads to a fine tuning 
problem: the relative values of p^ and pm niust be set 
to very high levels of accuracy (p0/pm)initiai ~ lO^^'^'^ in 
order to ensure p^/ Pm ^ 1 at precisely the present epoch. 
A more reasonable assumption might be if the energy 
density in the 0-field were comparable to that of radiation 
at very early times - say at the end of infiation [|6| . This 
might even be expected if the 0-field were to be an in- 
fiationary relic, its energy set by an equipartition ansatz. 
However for the 0-field to remain subdominant until re- 
cently its energy density must decrease rapidly at early 
times. Such behaviour clearly cannot arise for polynomial 
potentials V{<j)) oc (p^ , < p <, few, for which p^ will 
rapidly dominate the total density resulting in a colossal 
A-term today if p^ ^ pi-ad initially. Fortunately there do 
exist families of potentials for which the behaviour of p^ 
is more ffexible. To illustrate this, consider a minimally 
coupled scalar field rolling down the potential 



V{(i))^Vo{cosh\4>-\Y. 

V{(t>) has asymptotic forms: 

V{(l)) ~ Voe-P^^ for |A0| > 1 {4 
Vicj)) ~ Vb(A</))2P for |A0| < 1 



<0), 



(2) 



(3) 
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where Vq = Vo/2p. Scalar field models with the poten- 
tial V{(j)) oc e~P^'^ have the attractive property that the 
energy density in (jj tracks the the radiation/matter com- 
ponent as long as the value of cj) is large and negative, so 
that §: 
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(wb = 0, 1/3 respectively for dust, radiation). During 
later times the form of V{(f>) changes to a power law (H) 
resulting in rapid oscillations of (f) about = 0. The 
change in the form of the scalar fleld potential is accom- 
panied by an important change in the equation of state of 



the scalar field. As long as V{(j)) is described by (|^), the 
kinetic energy of the scalar remains larger than its po- 
tential energy ^0^ > V{(j)) and the scalar field equation 
of state mimicks background matter w^ ~ wb ■ However 
during the oscillatory phase {^(p^) can become smaller 
than {V{(j))), the virial theorem then gives the following 
expression for the mean equation of state |8] 
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The corresponding value of the scalar field density and 
expansion factor is given by 



,-3(l+ui^) 



acct", c= 1(1 + (w^)) ^. 



(7) 
(8) 



From (g), (0) & (§) we find that the mean equation of 
state, the scalar field density and the expansion rate of 
the universe depend sensitively upon the value of the 
parameter p in the potential (||). Three values of p should 
be singled out for particular attention since they give rise 
to cosmologically interesting solutions: 

1. p — 1: In this case the scalar field equation of state 
behaves like that of pressureless (cold) matter or 
dust (Wfj,) ~ 0, a scalar field potential with this 
value of p could therefore play the role of cold dark 
matter (CDM) in the universe. 

2. p = 1/2: This results in (w^) ~ —1/3 and p^ oc 
a^^. This choice of the parameter p leads to a 
'coasting' form for the scale factor at late times: 
a{t) oc t. A flat universe under the influence of the 
potential ¥{(/>) = Vo(coshA0 — 1)^/^ will therefore 
have exactly the same expansion properties as an 
open universe without being plagued by the 'omega 
problem' ! (See Q for related scenarios.) 

3. Smaller values p < 1/2 lead to {w^) < —1/3. From 
(0) we find that the scalar field density falls off 



slower than either radiation {pr 



matter {p„ 



,-3 



) or cold 



). The scalar field therefore 



dominates the mass density in the universe at late 
times leading to accelerated expansion according 
to a ex i^ c = 2/3(1 + (w^)) > 1. The epoch of 
scalar field dominance commences at the cosmolog- 
ical redshift 1 + z^ ^ (J^^/rj,^)-!/^^"'-/-). 

We therefore find that a scalar field with the potential 
(pi) might serve as a good candidate for quintessence. Fig- 
ure 1 confirms this by showing the density parameter fi^ 
as a function of the cosmological scale factor. We find 
that the ratio Q^/fls remains approximately constant 
during the prolonged epoch of radiation/matter domina- 
tion as the (/)-field tracks the dominant radiation/matter 



il 



1.00 



0.80 



0.60 



0.40 



0.20 



0.00 
10 



^ . 


/ ^ \ 


/ \ 


■^ \ - 


\ 


\ ' \l 


1 \ 1 


— "» "■■■ 1 H 


a \ 1 \ 


---a \' 


"cDM y ( 


A A 


/ \ / 


/ ^ / 


' \ / 


/ \ / - 


'' ■■■ / 


/~~~^ \ / - 


/ 1: J 


^ ■■■-.. 







10 10 10 10 10 10" 

scale factor a 



10" 



FIG. 1. The evolution of the dimensionless density pa- 
rameter for the quintessence field fi^ (solid line) is shown for 
the potential (Q) with p = 0.2. Matter (dashed line) and ra- 
diation densities (dotted line) are also shown. For z 3> 1 the 
quintessence field satisfies the tracker solution (|5|) and con- 
tributes a fixed fraction to the total density of background 
matter/radiation. At later times z ^ 2 scalar field oscilla- 
tions commence and the density in the quintessence com- 
ponent rapidly dominates the mass density of the universe 
leading to fi^ ^ 0.7 today. In this figure the parameters of 
the quintessence model used are: Vo = 8 x 10~^m%iMpc~^ , 
X = SOmT,} . 



component. [Q^/fls < 0.2 is necessary in order to sat- 
isfy nucleosynthesis constraints.) At the end of the mat- 
ter dominated epoch Q^/fls begins to grow as the scalar 
field equation of state turns negative in response to rapid 
oscillations of (f). As in earlier tracker quintessence mod- 
els Q present-day values fi^ ~ 0.7, ^m — 0.3, consistent 
with supernovae |^,^ and other observations |^, can be 
obtained for a large class of initial conditions. (General 
potentials leading to oscillatory quintessence must satisfy 
{V — (pV) > where angular brackets denote the time 
average over a single oscillation |1^. We assume that 
the quintessence field couples very weakly to other mat- 
ter fields so that rapid oscillations of cj) do not result in 
particle production of the kind associated with 'preheat- 
ing'.) Our potential belongs to the general category of 
exponential potentials which are frequently encountered 
in field theory fT^] , condensed matter physics |12] and as 
solutions to the non-perturbative renormalization group 
equations |13| . 

In figure 2 we compare the redshift dependence of 
the luminosity distance for a specific realisation of our 
quintessence model with cIl obtained from supernovae 
observations. 
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FIG. 2. The luminosity distance vs. redshift for the model 
shown in fig. 1 (solid line). The dashed line is the standard 
CDM model (shown for comparison) and the horizontal line 
corresponds to the fiducial empty Milne universe Qm -^ 0. 
The filled circles show supernovae data from High-Z Super- 
nova Search Team [^ and the opaque circles show data from 
Supernova Cosmology Project M. The low-z supernovae are 
from the Calan-Tololo sample. 

We would also like to draw attention to the possibility 
of a unified picture of quintessence and cold dark mat- 
ter in which both components are described by a pair of 
scalar fields evolving under the action of the potential (g) 
but with different values of the exponent p: 

y(0» = V^{coshX^c^-l)P*+V^,{coshX^^-l)P^ 

(9) 

p^ ~ I in the case of CDM and p^ < 0.5 in the case 
of quintessence. This approach (along to the lines sug- 
gested by |14)) ameliorates the 'coincidence problem' be- 
tween dark matter (CDM) and quintessence which arises 
in standard cosmology. It also significantly reduces the 
discrepancy between the present value of pm / Pr and that 
at the end of Inflation. Figure 3 shows a working exam- 
ple. 

It is interesting that the CDM particle in this scenario 
can be ultra- light, its mass [m? = V" ~ iJ^(t*)) is re- 
lated to the epoch t^ when ip begins to oscillate and its 
Compton wavelength Ac = m~^ can easily be of order a 
kilo-parsec or smaller. (In the cosmological model illus- 
trated in figure 3 the CDM field ■0 begins to oscillate at 
z « 10^ so that m~^ ~ 228 parsec.) Cold dark matter 
made up of a condensate of ultra-light particles would 
be frustrated in its attempts to cluster on scales smaller 
than Ac because of the uncertainty principle, the result- 
ing Frustrated Cold Dark Matter model (FCDM) might 
provide a natural explanation for two major difficulties 
faced by the standard CDM scenario, (i) The dearth of 



halo dwarf galaxies: the number of dwarf's in the local 
group is an order of magnitude smaller than predicted by 
N-body simulations of SCDM |1^]. (ii) The discrepency 
between observed shapes of galaxy rotation curves and 
simulated dark matter halos. Recent observations of low 
surface brightness (LSB) galaxies show them to possess 
rotation curves which indicate a constant mass density 
in the central core region. These observations are diffi- 
cult to accomodate within the SCDM model since high 
resolution N-body simulations of SCDM halos indicate a 
cuspy central density profile having the form p oc r 
in the core region |l6-20|. 
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FIG. 3. The evolution of the dimensionless density param- 
eter for the CDM field Q^ (dashed line) and quintessence field 
Q^ (thin solid line). Baryon (dash-dotted line) and radiation 
densities (thick solid line) are also shown. The parameters for 
<j> have been given in the caption of figure 1. The parameters 
for Tp are: p^ = 1, V^ = 3 x l(fm%iMpc~^ , \^ — 8ml 
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It is both interesting and revealing that a physical 
mechanism suppressing small scale clustering arises in 
the FCDM model at the purely classical level. To demon- 
strate this we note that once -0 begins oscillating V{^) 
acquires the form T/(V') c^ Vol^X'^ij^ + ^A'^V'I. Us- 
ing the relationship ^ (i;^) = Vo[X'^{i^^) + ^(V''*)] 
and assuming that the gradient energy is subdominant 
we obtain the following expressions for the mean pres- 
sure p^ and the mean density p^: p^ — ■^A'*('0^), 
p^l, ~ VoA^('(/'^). The motion of ip is driven mainly by 
the quadratic term X'^ip'^ (which also provides the dom- 
inant contribution to the energy density). As a result 
ip = y/VoX'^{ipQ — ip"^), which can be used to establish 
(■04) — |(^2^2_ Substituting the resulting expression for 
the speed of sound v1 = dp^/dp^ = Pip/SVo into the 
Jeans length Aj ~ ^JvlJI-nGp^, we get 



A ,m,pi. 



.p, ^ ^( — ). 



2^0 V2 rn 



(10) 



We therefore find that the Jeans length is larger than 
the Compton wavelength (Aj > m~^) if A > rripj [ p2[ . 
For the FCDM model illustrated in figure 3, one finds 
Aj ~ f kpc. A lighter particle will possess a larger Jeans 
length while Xj is smaller for a V'-field which began oscil- 
lating much before t^q, a very massive field will resemble 
standard CDM. By inhibiting gravitational clustering on 
scales smaller than Xj ~ kpc, FCDM is expected to give 
rise to galaxy halos which are less centrally concentrated 
leading to better agreement between theory and observa- 
tions. 

Finally we note that quintessence-type potentials could 
also arise in particle physics models which invoke the 
Peccei-Quinn mechanism to solve the strong CP prob- 
lem in QCD. Consider for instance the following simple 
modification to the symmetry breaking potential respon- 
sible for the axion 



y(</>) = xM' 



f 



r+m'ni 



(11) 



The second term in ( [11[ ) when expanded about 9 = 
a,rg{(j)) = acquires the form m?p{9'^ /2Y. Accord- 
ingly rapid oscillations of the field about ~ Q now 
give rise to an equation of state described by (|6|), re- 
sulting in {wg) < —1/3 for p < 1/2. (In the standard 
scenario p = 1, (we) — 0, m ^ rua is the axion mass 
and / -^ fpQ is the Peccei-Quinn symmetry breaking 
scale.) As a result an axion-like scalar with p < 1/2 will 
be a candidate for quintessence since its energy density 
will diminish more slowly than that of either matter or 
radiation, leading to the dominance of the 0-field at late 
times and the accelerated expansion of the universe. 

It should be noted that motion under the action of po- 
tentials (||) & (|ll|) is well defined even though for p < 1/2, 
V'{(l)) is weakly singular at = 0. One can make V{(f)) 
mathematically more appealing by the field redefinition 
(/) — + (0^ -I- (piy^^, 4>c -^ 0, this will not affect our results 
in any significant way. 

It is also worth mentioning that V" < during the 
oscillatory stage if p < 1/2. This is likely to affect very 
long wavelength fluctuations in for the potential (g) 
since the scalar field begins oscillating fairly recently in 
this case. 

Finally we would like to point out an important dis- 
tinction which exists between the quintessence models 
suggested by us and those of Q. In our models the 
quintessence field (p ov oscillates about a small value 
at late times (formally (/), 6* ^ as i — + oo). Additionally 
the potential V{4)) is not constraind to be flat since the 
field does not have to slow-roll in order to give rise to a 
negative equation of state. As a result quantum correc- 
tions which might be significant in quintessence models in 
which the scalar field rolls down a flat potential to large 
values (j> > mpi |Q can safely be ignored in models of 
the kind discussed in the present paper. 
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